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Abstract. We present a method to compute the holonomic extension of a D- 
module from a Zariski open set in affine space to the whole space. A particular 
application is the localization of coherent D-modules which are holonomic on 
the complement of an affine variety. 



Throughout this article let R — K[x±, . . . , x n ] be the ring of polynomials in n 
variables over the field K of characteristic zero and D = R{d\, . . . ,d n ) the n-th 
Weyl algebra over K. Here, di — Suppose / G R, let J be a left ideal of D 
and set M = D/J. Let X = K n , Y*= V(/) C X, U = X \ Y and j : U ^ X the 
natural inclusion. 



In the landmark paper K2] M. Kashiwara proved that J j M is holonomic 
provided that M is holonomic on U, but not necessarily on all of X. In this note 
we make this important result algorithmic in the sense that we provide an algorithm 
that computes the module structure of J j _1 M over D. 



We remark that U is affine and X = K n , and so J-j M = Rf ®r M as D- 

module where Rf = -R[/ -1 ] = T(Ox, U). An algorithm to compute Rf ®ij M as a 
D-module has been given in [ Oa3| un der the assumption that M is holonomic on 



K n and /-torsion free, and in ]Q-TS| under the assumption that M is holonomic 
on K n . The advantage of our new algorithm, besides removing the requirement 
of holonomicity on V(/), is that the natural map M — > Rf <E)r M can be traced, 



which is an important property for other computations as local cohomology ([Oa3|, 



m and dc Rham cohomology ([O-Tl], [W2|) 



The algorithm is illustrated in a section devoted to computational examples. 

1. Description of the Algorithm 

We shall start with a motivating construction which will also provide a skeleton 
of the proof of the correctness of our algorithm. 

Assume that K is in fact the field C of complex numbers. Let s, v be two new 
variables and let D v = D(v,d v ). 

Notation 1. Throughout we will denote df/dxi by /j. If P G D v is an operator 
interpreted as a polynomial in 8%, . . . ,d n , we will simply write P(d x ). Similarly we 
will abbreviate P(di - v 2 fid v , ... ,d n - v 2 f n d v ) by P(d x - v 2 f x d v ). 

Consider the ring D v and the C-vectorspace Rf[s] <E> f~ s - Here, / _s should be 
thought of as a symbol that behaves like the complex valued function /~ s under 
differentiation and the tensor is over C. We shall turn this vectorspace into a 
module over D v by means of the action • defined as follows. We require that 
all Xi act as multiplication on the left factor, while <9, acts by the "product rule" 
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di • {^P 1 ® f- s ) = dii^pl) ® /— + g /-«. Thus, we only need to 

define the action of v and d v . We set 



f m J p 

We observe that v is simply shifting s up by one (in particular, the extra / in the 
denominator comes from shifting the exponent in / _s ) and d v is shifting s down 
by one, together with the "differential" 2 — s of the shift operator (cf. the action 
defined in Q). 

Armed with this action on Rf[s] ® f~ s we can define an action of D v on Dj J ®r 
{Rf[s\ (g) /~ s ) as follows, v and d v commute with the first factor in the tensor 
product, Xi acts by left multiplication on the first factor and di acts by the product 
rule. 

We detect a certain set of trivial relations for this action. Namely, let us ask 
which operators on D v are candidates for annihilating the element 1 (g) 1 <g> f~ s , 
the bar denoting cosets in D / J . Independently of J, 1 — vf is such an annihilator. 
Assume that J is generated by {Pj{d x )}i where we think of Pj as a polynomial in 
d x with left coefficients in R. Let us define a ring map <f> from D to D v by mapping 
P G D to <f>(Pj(d x )) = Qj(d x ) := Pj(d x - v 2 f x d v ) G D v (recall /< = $(/)). It is 
easy to see that this assignment preserves commutators and hence is an actual ring 
map. We remark that (j) represents the map U V(l — vf) C C n+1 given by 
x — > (1/ f(x),x) on the level of differential operators, while the inclusion </>(-D) 
Z?„ corresponds to V(l — vf) C" +1 . 

We claim that Qj(d x ) annihilates 10 1® / _s . To prove this, we observe that 
{di - v 2 fjd v ) • (P® 1 ® equals ® 1 ® /" s for all P E D. Thus Q^) • (T® 
1 ig) / _s ) = Pi{d x ) ® 1 <8 / _s = 0. We point out that — + 1 acts as multiplication 
by s and that this construction is (formally) independent of the base field. 

Our main statement is the following 

Theorem 2. Let M = D/J. Consider the left ideal D v ■ {1 — xv,4>{J)} and the 
right ideal {d v } ■ D v in D v . The quotient D V /{D V ■ {1 — xv, 4>{J)} + {d v } ■ D v ) is 
isomorphic to M <g> Rf as D-module. Moreover, the natural map M — > M ®r Rj 
sends 1 G M to 1 (g) -p. 

It is worth noticing that this is in fact a left D-module due to the fact that D 
and d v commute. 

In the general context in which the theorem is stated, the D-module D V /{D V ■ 
{1 — xv, <j){J)} + {d v } ■ D v ) will not be a finitely generated D-module. If however M 
is holonomic on[/ = K n \V{f), then M® R R f = D v /(D v -{l-xv, cf>{J)}+{d v }-D v ) 
will be holonomic, and in particular finitely generated. In order to compute the 
structure of the localized module, we observe that it is the integration of the left 
D„-module D v /D v -{l — xv, <fi{J)} with respect to d v (cf. Oa3, Theorem 5.7], |Q-T1 
Section 6]) 



Methods for the algorithmic computation of 0-th integration worked out in ]Oa3 



Theorem 5.7] (i.e., the restriction to v — of the Fourier transform of the module,) 
together with our theorem give the following algorithm. Let N denote D v /({d v } ■ 
D v +D v -{l-fv,<f>{J)}). 
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Algorithm 3. Input: / G R; {Pi, . . . , P r } C D generating J, D/J holonomic on 
U = K n \V(f). 

Output: k S N; {Qx, ■ ■ ■ , Qt} Q D such that D/D ■ {Q%, . . . , Qt} is isomorphic to 

D/J (g>fl i?/ generated by 1 <g> pW- 

Begin. 

1. For Pi,i — 1, . . . , r compute </>(Pi) £ defined by replacing d x by d x — 

2. Compute the 6-function b(s) for integration of D v /D v ■ {1 — fv,<fi(J)} with 
respect to 9^. That is, find K[vd v ] (~l in u) (D lI • {1 — (j)(J)}) where w is the 
weight assigning 1 to v, —1 to d v and to all other variables. Replace vd v by 
-s - 1. 

3. Let k be the largest non-negative integer root of b(s). If there is no such root, 
then output I = D, else continue. The integral N is generated by the coset 
of v k in N. 

4. Compute the annihilator / over D of v k in N. 

5. D/J ®r Rf is generated by 1 ® /~( fc + 2 ) and isomorphic to D/I. 

6. Return k and /. 
End. 



An algorithm to perform step 2 of the algorithm is given in |Oa2, Algorithm 4.5]. 
The steps 3 and 4 are nothing but an algorithm to get the 0-th integral of D- 
modules. Here is a more precise description of steps 3 and 4. 

1. Let G = {g\, . . . , g m } be a Grobner basis of D v ■ {1 — vf, 4>(J)} with respect 
to w (see, e.g., [ BSTj , section 1.1]). 

2. Let Gt be the set 

{ normalForm^gj, {d v } ■ D v ) \j = 1, . . . , m, < i < k — ordu,^)}. 

3. Regard D ■ {Gk} as a left submodule of the free module 

D v° + D v 1 + --- + D v k = D k+1 

and find generators of D ■ {Gk} flfl v k . The last intersection can be computed 
by using an order to eliminate v°, . . . , v k ~ 1 , i.e, by using an order y such that 
CUV 1 y akv k for all i = 0, . . . , k — 1 and a,i, at ^ in D. 
Here, we put 

ord w (g) := max ju; • (a,/9) | g = ^a af3 v a d^ 

following the notation of ]SST ]. normalForm(<7, {d v } ■ D v ) means taking the normal 
form of g with respect to the right ideal {d v } ■ D v . For example, 

normalForm(> 2 <9 2 , {d v } ■ D v ) = normalForm(<9 2 v 2 - Ad v v + 2, {d v } ■ D v ) = 2. 

2. Proof of correctness of the algorithm 

Let us first provide the 
Proof of the Theorem: Given an operator P 6 D v we shall define its normal form 
in D v . The goal is a presentation of P where di — v 2 fid v take the position of di. In 
other words, we are aiming for a sum of the form P = J2 dv vb Pab(x)qab{d x —v 2 f x d v ) 
where p a b and q a b are polynomials in n variables over K. To this end, write first 
P as Y l d%v b pab(x)q a b{d x ). The operator P 1 := J2 d v vb Pab(x)q a b(d x - v 2 f x d v ) has 
the property that P — P 1 will have lower degree in d x than P. If P — P 1 = 0, quit. 
Otherwise write P - P 1 as £ fljAiiC*)^^) and set P 2 := S^^^x)^^^ - 
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v 2 f x d v ). Repeat this procedure until we arrive at P l = 0. Then P = P 1 + ■ ■ ■ + P l 
and this sum is the desired normal form in D v . We shall write P for the normal 
form of P in D v . 

The normal form in D v of an operator induces a normal form in D v / {d v } ■ D v by 
removing all terms in {d v } ■ D v from the normal form of P in D v . We denote the 
normal form of P G D v in D v / {d v }-D v by P. Of course, P+{d v }-D v = P+{d v }-D v 
as cosets. We notice that both normal forms are unique. As an example, consider 
the normal form of <9j in D v /{d v } ■ D v : P = di, P 1 = di — v 2 fid v , thus P 2 = 

v 2 fid v = d v v 2 f l - 2fiV, P 3 = and hence P = (d t - v 2 fid v ) - 2f t v + {d v } ■ D v . 

As a D-module, the quotient N = D v /({d v } ■ D v + D v ■ {1 — vf,<f>(J)}) is 
spanned by the powers of v, which are in normal form. Let us try to understand 
the relations among these powers. To this end note that (1 — fv)(pi — v 2 fid v ) — 
{d l - v 2 fid v + fiv)(l - fv) which implies that (1 - fv)P = in N for all P G D v . 

This proves that f(vP) = P = P in N for all P <G D v and we conclude that every 
element in N is divisible by /. 

Moreover, the division by / is unique since if fQ = fQ' in N then vf(Q — Q') = 
in N and that implies Q = Q' in N. Conversely, for all P e N there is a 
unique Q e N for which vQ = P because then necessarily Q = fP. So / is an 
invertible operator on N and the inverse is given by v acting on the normal form 
representative of a coset. 

We shall now make use of the action we defined in the previous section. 

The map D v -> D/J (R f [s] ® /^ s ) sending 1 to T ® 1 (g) f~ s defined via the 
action • is surjective. Consider the combined map D v —> D/ J ®r {Rf[s\ ® f~ s ) — > 
(D/J ® R (R f [a] ® r s )) /{a - 2) (D/J ® R (R f [s] ® f- s )). 

This is a surjective map of left D- modules and the kernel contains {d v } ■ D v + 
D v • {1 — vf, 4>(J)}. If an operator in D v is in the kernel, we can assume it to be in 
normal form in D v /{d v } ■ D v since {d v } ■ D v is in the kernel. Any such operator may 
be written as v k P(d x — v 2 f x d v ) where P(d x ) E D. We conclude that P must already 
be in the kernel since v k P(d x - v 2 f x d v ) • (T <g> 1 ® /" s ) = v k • (P(d x ) 01® / ~ s = 
~P(d^®f- k ®r s . Thus P(d x ~v 2 f x d v )*(T®l®r s ) e (s-2)D/J® R (R f [s]®f- s ) 
and therefore fiP(d x ) G D[s] ■ (J, (s — 2)) for some j G N. Since neither P nor 
J does contain s, setting s = 2 we see that P(d x ) is in the /-saturation of J, 

pp(di) g j. 

PP(d x ) G J implies that PP(d x - v 2 f x d v ) G 0(J). Modulo D„ ■ {1 - vf} 
this is equivalent to P(d x — v 2 f x d v ) G 0(J) since P(d x — v 2 f x d v ) is in normal 
form. This simply means that P is in the kernel of D v — > D/J (& R (Rf[s\ ® 
/" s ) -> (^/J®i? (i?/[s]®/- s ))/( S -2) (D/J® fl (i?/[s]®/" 5 )) if and only if it 
is in {94-C + C-{l~w/,0(J)}. 

Finally we observe that 

(D/J ® R (R f [s] ® /"*)) /(s - 2) (D/J ® R (R f [s] ® /- s )) 

is isomorphic to D/ J® R Rf by means of the map P® ® f~ s ^P®^^. □ 

Remark 4. The theorem and its proof generalize nearly verbatim to the situation 
where M is any finitely generated D-module. 

Now let us consider the situation in which M = D/J is holonomic on U. Then 
D v /({d v } ■ D v + D v ■ {1 — fv, (j)( J)}) is a finite D-module as it is isomorphic to the 
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module M®nRf which is holonomic by theorem 1.3 of JK2| | . We can find a generator 
by computing the 6-function for integration along d v for D v /D v ■ {1 — vf, 4>{J)} as in 
step 2 of algorithm [| If k is the largest root of the 6-function then v k is a generator 



for D v /({d v } ■ D v + D v ■ {1 - fv,<j){J)}) = M ® R f (compare |0-T2|, algorithm 
5.4). Thus in order to represent the localization M Cg).R Rf as a quotient of D all 
one needs to do is to find the annihilator of v k over D. This shows the correctness 
of our algorithm. 

Remark 5. Again, the algorithm generalizes to the non-cyclic situation. Let M — 
D m /J, D' m = (B™Dej. The modifications are as follows. Compute m separate 
^-functions bj(s) to the integration of D v (ej + (j)(J)) / (D v ■ {4>{J), (l-vf)ej}) along 
d v . M®rRj is generated by the cosets of v kj ej in N where kj is the largest integer 
root of bj(s). 



3. Categorical explanation of the algorithm 

Let us now give a more categorical explanation of the validity of our algorithm. 
Decompose j : K n \ V(f) — ► K n as j = p o l o (j); 

W K n+1 

[' 

K n \ V(f) K n 

where : K n \ V(f) — > W — {(v, x) \ v ■ f(x) = 1} is defined by x i-> 

L '. W e — > K n+1 is the closed embedding, and p : K n+1 — > K n is the natural 

projection. Then we have 



M® R R f = / r L M = r l M 

J j J p J L J (j) 



by the chain rule of integration functors (see, e.g., |Bl p. 251, 6.4 Proposition], [HT 



1.5.1]). Note that the chain rule holds in the derived category in general, but, in 
this case, j. is an exact functor since T> } = T>x\\/ f] is flat over T>u and j is an 
affine morphism. J t is also an exact functor. Hence, we have only to compute 
0-th integrals in each step of the computation of the integral functors. /, j~ x M is 
obtained by the coordinate transformation represented by our ring map 0; consider 
the map 

Vi =Xi, (i = 1, ... ,n), y n+ i = l/f(x). 

Then, d Xi = d Vi — -j%dy t = d Vi — fiV 2 d Vi modulo 1 — fv = 0, which commute with 
each other and define our ring map <j>. J is nothing but Kashiwara equivalence cor- 
responding to 4>{D) ^ D v . It follows that f l Lj ~ 1 M = D V /(D V ■ {(p(J), 1 - fv}), 
compare also Proposition A.l (p. 596) of |Oal|| . Integration under the projec- 
tion p corresponds then to our last step accomplished by forming the quotient 

(.U,./ R!-(,U,./ 

Since j~ x M is holonomic, so is J t j,j M and hence the 6-function for inte- 
gration is nonzero, thus guaranteeing termination of the search for generators of 
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4. Examples 

Example 6. For our first example we take n = 1, x = Xi, K = C, J = xd x + A, A 6 

C,f(x)=x. 

In this scenario, we have to compute the integral of the module D v / (1— xv, x(d x — 
v 2 d v )+X) along d v . As D v -{1— xv, x(d x —v 2 d v )+\} = D v -{d x x—d v v+l+X, 1— xv}, 
one sees that the b- function is s(s + 1 + A). Thus the largest integer root is either 
or — 1 — A, depending on whether A is a negative integer or not. So M ®r R x is 
generated by 1 (g) a;~(~ 1 ~ A ) _2 in the former and 1 ® x~ 2 in the latter case. 

If for example A = —7 we compute a 6-function of s(s — 6) indicating that 
M <S)r R x is generated by 1 <S> a; -8 - Since in this case M = R generated by x 7 , we 
conclude that M <g>fl Rf is in fact generated by x~ 8 ■ x 7 — 1/x, as it should. 

If on the other hand A = 1/2 then b(s) — s(s + 3/2) and hence the largest integer 
root is 0. Thus, M is generated by the germ of x~ x and M ®rR x is generated by 
v° corresponding to l/x 2+x . M is already isomorphic to M ®r R x . 



Example 7. In this example we consider the left ideal J generated by 

d x (x 2 -y 3 ), 
d y (x 2 -y 3 ). 

These are annihilators of the function 1//, where / = x 2 — y 3 , but they do not 
generate the annihilating ideal (see, e.g., ]Oa2| ). The left ideal J is not holonomic 
since the characteristic variety of J is V(x 3 —y 2 )^JV(£ lX ,S l y), whose first component 
has dimension 3 in C 4 = {(x, y; £ x , £, y )}- As to an algorithmic method to get the 
characteristic variety, see J0a0 |. We give now the output of a computer session 
using the computer algebra system Kan/sml (Q) interspersed with comments. 
We remark that in this case M restricted to U is an Oy-coherent free module of 
rank one where U = C 2 \ V(/). 
/ff 

[ ((x"2-y"3)*(Dx - v~2*2*x*Dt) + 2*x) 
((x"2-y"3)*(Dy + v~2*3*y~2*Dt) -3*y~2) 
(v*(x"2-y"3)-l) 

] 

def 

(this is <j){J)) 



sml>ff [(v)] intbfm : : 

[ $216*s-4+1296*s-3+2586*s-2+1716*s$ ] 
[[6,1] , [s+2,1] , [6*s+ll,l] , [6*s+13,l] , [s,l]] 
(these are the factors of the 6-function) 
sml>ff [(v)] -2 1 intall_s ; 
Completed, 
(computing the integration) 
0-th cohomology: [ 1 , 
[ -3*x*Dx-2*y*Dy-18 , 
3*y~2*Dx+2*x*Dy , 
-2*y"3*Dy+2*x"2*Dy-18*y"2 ] ] 
-1-th cohomology: [ , [ ] ] 
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The integration J. j 1 M is not Ox -coherent, although of course it is still coherent 
over the sheaf of differential operators T>x on X = C 2 . Since localization is an 
exact functor, the first cohomology group (corresponding to the first higher order 
integration) was known to be zero. The 0-th cohomology above coincides with the 
annihilating ideal of the function /~ 3 . 



Example 8. Let n = 3, and consider the ideal J generated by the system 

(x 3 ~y 2 z 2 fd x + 3x 2 , 
J y 



(x 3 -y 2 z 2 ) 2 d y - 2yz 2 , 



(x 3 -y 2 z 2 yd z - 2y'z. 

These operators are annihilators of the exponential function e 1 ^ where f(x, y, z) — 
x 3 — y 2 z 2 . The characteristic variety of M = D/J has six components, defined by 
the prime ideals (y, x), (z, x), {£, Xl £ y , £z), {£,y, z, x), {£, z ,y, x) and the ideal generated 

by 

y£y - z£ z , 2x£ x + 3z&, 8z£ + 27^&, 

8yg + 27^1 -Az 2 e x + 9z£ 2 , -4yz£ + 9atf w &, 

-4y 2 £ 2 +9^ 2 , 2z 3 Uz+3x 2 e y , 2yz 2 £ x + 3x 2 £ y , 

2y 2 z£ x + 3x 2 £ z , -yz 3 i z + x 3 £ y , -y 2 z 2 + x 3 , 

x 3 e y -z^i 

All but the first two are of dimension three. This implies that the non-holonomic 
locus of M is contained in the hypersurface x — 0. Here, we used [Asir] to obtain 
the primary ideal decomposition. 

Hence we may apply our algorithm to compute M ®r R x . We remark that 
contrary to the previous example in this case j~ l M is holonomic but not coherent 
as Cc/-module on U = X \ V(x). Using Kan/sml again one obtains the eight 
operators 

-3yd y + 3zd z , -2xyz 2 d x - 3x 3 d y - 4yz 2 , -2xy 2 zd x - 3x 3 d z - 4y 2 z, 

6xz 3 d x d z + 9x 3 d 2 + 6xz 2 d x + 6yz 2 d y + 6z 3 d z + 12z 2 , 

-6y 2 z 3 d z + Ax A d x + I2x 3 zd z + 8x 3 + 12, 

6yz 4 d 2 - 4x 4 d x d y - 12x 3 zd y d z + I8yz 3 d z - 8x 3 d y - 12d y , 

8x b d 2 x + 2Ax A zd x d z + I8x 3 z 2 d 2 z + Ux A d x + I02x 3 zd z + 80x 3 + 2Axd x + 48, 

-6z 5 d 3 + Ax 4 d x dl + 12x 3 zd 2 d z - 36z 4 <9 2 + 8x 3 d 2 - 36z 3 d z + 12d 2 

* y y * y ~ y 

which annihilate the function x^ 2 e x ^ . The characteristic variety of this holonomic 
left ideal of D has the same last four components as the characteristic variety of 
M while the first two components (of dimension 4) are replaced by (^ z ,(, y ,x) and 
(z, y, x) (of dimension 3). 
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